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9. The Secondary Field considered as the Superposition 
of an Infinity of Waves originating after 1, 2,3.... 
Reflections. 


i is Part I. of this paper the diffraction round a sphere 

of waves originating from a point source was 
considered starting with the harmonic series (1). There- 
upon a transformation of this series was affected with. 
the aid of a contour integral, resulting m a series which 
converges rather rapidly in the shadow region, while the 
convergence in the “lit ” part of space is not appropriate 
for numerical evaluation. As in many physical applica- 
tions a knowledge of the wave-field in the “ lit ” part of 
space is of importance, we now proceed to an entirely 
different attack on the problem; while again starting 

* Part I. of the present paper appeared in the Phil. Mag. xxiv. 


p. 141 (July 1937). 
+ Communicated by the Authors. 
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with the rigorous series (1), the fundamental idea used 
here is to split up each of the terms of (1) into a geo- 
metrical series, after which the corresponding terms of 
these geometrical series are combined. This splitting- 
up process is based on the physical conception of radial 
waves being reflected several times between the surface 
and the centre of the sphere. Thereupon each of these 
new series of radial waves is written in the form of a 
multiple exponential integral. It next appears that 
the position of the saddle points of these integrals 
determines the geometrical points at the surface of 
the sphere, where a ray joining the point source and 


{a) infinite slab, (6) sphere, showing the different elementary waves. 


the observer is refracted and reflected. Thus our saddle 
point treatment of wave functions appears to be identical 
with the application of Fermat’s principle. The second 
order approximation of these integrals with the method 
of stationary phase yields the amplitude of these rays 
whose geometrical form was given by the above-mentioned 
saddle points. These rays correspond exactly with the 
classical geometrical ray treatment of the rainbow which 
originates from Descartes. 

In order to illustrate the fundamental idea of the 
splitting up of the wave function in successively refracted 
and reflected waves we treat in this paragraph an 
analogous simple problem. 
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Consider the problem with two plane boundaries repre- 
sented by fig. 12a, where an electrical vertical dipole is 
situated at z=b, p=Oin the medium 1 with propagation 
constant k. Medium 2, with propagation constant ka, 
consists of an infinite slab with plane boundaries at z=a and 
z=0. Below z=0 we assume a third medium of infinite 
conductivity (k,= ), so that all waves penetrating into 
the second medium are totally reflected at the boundary 
surface z=0. Later on we shall see that in the spherical 
analogy (fig. 126) we again find a total reflexion, but there 
against the centre of the sphere. This plane problem is 
easily solved with the well-known integral method of 
Sommerfeld, and leads, for a primary field e*'®/(7k,R), to a 
secondary field in the first medium, given by 


i — z 
Tec = Æ y .RA)-Ja(Àp). e— Vx k(e+b—20) 


7 
(60) 
where 
V¥—k,? k? Ne E plete 2 Tgh Rh Ea) 
k 
Pp arts ars ky? - ae “9 (61) 
be + VANE tg .Tgh( B—k,?a) 


and p is the horizontal distance from point source to 
receiver. 

If, with Sommerfeld and Weyl, (60) is considered as 
the superposition of plane waves, A has the meaning of 
k, sin $’, where 3’ is the angle a particular plane wave 
makes with the vertical. The fundamental idea of this 
paragraph is to write (61) as the geometrical series 


R,=RQ)+ F (1+RQ)} 
-RA)E{1—R(A)} .e- 28 Nav ie (62) 


so that, when inserted into (60), the secondary field in the 
first medium becomes the sum 
df. =02,+0,°+0,/°+0,'+ 0,°+ er . (63) 
3H2 
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where all the 0% mean outgoing waves in the first medium, 
represented by * 


ee | Add ORAA) Jo(Ap) 67 VEFRe+ 20), 


ik, Jo Bk? 
(64a) 
. 1 f? Ada 
o= j are {14+B(a)}.R(A)*™.{1—R(a)} 


J (Ap) eV bret b—24)—- 21K +a V A-k (K>0), (645) 
and where R (A) has the usual form 
7 ie aie a A 
kv Bh +-k eV Pk? 
meaning physically the Fresnel reflexion coefficient for 


the waves incident at an angle 9’, cf. (21). The 


physical meaning of (645) becomes at once clear if in 
the integral we put 


A=k, sin 3’=k, sin 3”, 


(65) 


so that we obtain 


zjį2— io 

k= f” sin s d9".(1+-RQ)}.ROVE.(1-RQ)} 
0 

-J (Ap) etka cos $'(z-+b—2a)+ 2ike(K+1)a cos 9” | (66) 


The form of (66) is easily interpreted as an outgoing 
wave in the first medium (see fig. 12 a), originating from 
the primary incident wave I’, in the first medium, which 
first is refracted as the ingoing wave I‘ in the second 
medium, this being afterwards totally reflected at z—0, 
thus becoming Of; this outgoing wave, when reaching 
z=a, is split up in the refracted wave O¢ in the first 
medium and the reflected ingoing wave Ii in the second 
medium. This process is repeated till finally we obtain 
Oz as the outgoing wave in the first medium, resulting 
from waves in the second medium which have been 
reflected up and down in the second medium with K 
reflexions at z=a and K+1 at z=0. It will be clear 
that the total secondary field Heo in accordance with 
(63), is obtained by superposing all the outgoing waves 
characterized by K=—1,0, +1, +2.... ete. in the 


* For the symmetry of the notation we give the incident and first 
reflected wave the suffix —1. 
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first medium. The notations are chosen such that 
O and I relate resp. to the outgoing and ingoing waves, 
whereas the index e or i denotes the first or second 
medium. 

The above physical consideration now makes the form 
of (66) quite obvious, because the integrand contains 
the factors 


{1-4 R(A)}.R(A)*.{1—R()} 


2ike(K+i1 ” 
and erika K + i)a cos $7. 


the first of which represents the attenuation 1+R(A) at 
the first refraction, the attenuations R(A)* due to the 
K reflexions at z=a, and the attenuation 1—R(A) at the 
final refraction into the first medium. 

The second factor represents the phase retardation 
and, for complex k,, the attenuation due to the fact 
that the wave has travelled K+1 times up and down 
in the second medium, the form of this factor being 
quite similar to that obtained in the derivation of the 
phase retardation in Bragg’s well-known formula of 
crystal reflexion. 


10. The “ Spherical Reflexion Coefficients.” 


In the last paragraph when we split up the secondary 
field into an infinity of components we were naturally 
lead to introduce the ordinary Fresnel reflexion coefficient 
which belongs to the reflexion of a plane wave against 
a plane boundary. Before considering in more detail 
the analogous splitting up of the waves in our spherical 
problem we shall have to investigate first the extension 
of the simple Fresnel reflexion coefficient for the reflexion 
of an incident spherical wave 


C@ (kyr).P,, (cos 3) 
against the concentric sphere r=a separating the first 


medium (k,) from the second medium (k3). 
We can write the reflected wave in the first medium as 


1) 
Ry, -£(k,a). rane .P (cos3),. . (67a) 


and the refracted wave in the second medium as 


2 
Ryo. t2 (ka). Fae .P,,(cos3). . . (676) 
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These forms have been chosen such that, e. g., R,, means 
the ratio of the amplitudes of the reflected and incident 
waves at r=a, and similarly for R,.. 


The boundary conditions, A (rI) and k?IT continuous 


at r=a, now fully determine the “ spherical reflexion 
and refraction coefficients”? R,, and R,,, which in our 
case, where we assume the electric field to be in the plane 
of incidence (electrical dipole}, become 


Rau = 


[iaoee] + nee], 
Ee T log EZA vet] ` — z Ai log eTA = 
(68a) 
Re =gh +Rn). - aaa a ee ee (688) 


It is of interest to contrast the more physical factor 
R,,, which contains ¢ functions only, with the formally 
used R,, of (2), which contains both ¢ and # functions. 


ll. The Splitting-up of the Secondary Field into an 
Infinity of Spherical Waves originating after 1, 2, 
. Reflexions at the Centre and the Surface of the 

Sphere. 


After having derived the ideas of the spherical reflexion 
and refraction coefficients we are now able to interpret 
the spherical case as an analogue of the plane problem 
of §9. To this end we follow again the waves on their 
courses. According to (36) the primary field for r<b 
contains the standing waves %,(k,r)P,(cos $). As 


tn(kyr) = 301 kır) ee (kyr), 


we can consider the primary field as containing both 
outgoing ¢@) and ingoing ¢ waves. The outgoing part 
will be presently left aside. The ingoing part is given by 


= p '(2n-+1). L2(k,b) LP (kir). Pa(cos 3). 


As CEN in § 10 each of the components, say 
tP (k,r)P,,(cos$), is first at r=—a split up into a reflected 
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wave in the first medium and a refracted wave in the 
second medium, with the reflexion and refraction 
coefficients R,, and R, resp., thus, after summation 
over n, giving rise to the waves Oĉ, in the first medium 
and I? in the second (see fig. 126). Each of the com- 
ponents of the latter, say 


2) Aly (kar) 

H(2m+1) Ryo. OE bE (ya). ? SE) .P, (cos 3), 
would, however, cause a singularity at r=0, due to the 
factor (2 (kyr), which can only be compensated by an 
outgoing wave 


12n +1). Ria. (k b)i (ka) bn (kar) P, (cos 3) 

25 ->12 Sn V1 BN is Olka) sAn 3 
having the same amplitude. This fact can be interpreted 
as a total reflexion of the ingoing wave Ii against an 
infinitely small sphere at the centre. This is, moreover, 
confirmed by the fact that the limit for a=0 of the 
spherical reflexion coefficient (68 a) becomes unity inde- 
pendently of the values of k, and &,. Hence it can be said 
that the incident wave in the sphere is totally reflected 
against its centre analogous to the reflexion against the 
plane z=0 of the plane problem of § 9. 

The summation over 7 of these totally reflected waves 
results in the outgoing wave O}, but still being in the 
second medium. Following this wave further, it will 
arrive at the inside of the boundary r=a, where it will 
be divided into an ingoing reflected wave Ii in the 
second medium (see fig. 12 b) and an outgoing refracted 
wave O* in the first medium, where it has to be noted 
that now this reflexion and refraction occur with the 
factors Rə, and R,,, derivable from R, and Ris by 
interchanging k; and kə This process repeats itself 
indefinitely in a fully analogous way as described in the 
plane problem of § 9 (compare fig. 12a and fig. 125); 
therefore we have chosen the notations for the different 
waves quite the same as in § 9. 

We thus obtain instead of (63) and (64) the following 
final formule : 

for the total field in the first medium 


lHa = E H0 , +0,°+0,°+0,°+0,+...(r>a), (69a) 
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and similarly for the total field in the second medium 
Moa = Ió +Oo'+1,'+ 0; +1,'+-0,'+-1;'+05'+ . . .(r<a), 


(69b) 
where 
bE (2n-+1)L2(byb)L (kyr )Pn (eos 3), (>b), 
i, = : (70a) 
bE (n+ EME DIP ErP, (008 3), (r<), 
(70b) 


O=} È (2n -+R aie (kbila). S264"). P, (cos 5) 
1 2 pani Iln 1 n 1 CD(k a) 7 
(70 c) 
Ok = 4 È (n+ DR RER. EEDA) 
A EO (kr) 
` Taa) ` LO (ka) 7 


Tg =4 Š (204 DR RE. Pb Ea) 


P,,(cos$), (K 20), (70d) 


£0 (kya) E CO (ker) 
l E "Tk ya) 


Ok =} X (2n+1)R RE. CD(k b) (ka) 
n=0 


kee) E, DE) 
USP (ka) OP (kaa) 
In deriving these formule we started, according to (70 b), 
with the ingoing part only of the primary field. It can, 
however, easily be shown that, in order to satisfy the 
boundary conditions, the other, outgoing, part of the 
primary field must be accompanied by a secondary field 
which just cancels it, so that our solution (70) gives 
the total field, satisfying the boundary conditions and 
having as only singularity the point source at Q. 

It may be finally remarked that the physical arguments 
used in this section are equivalent to the development of 
the denominator of the general R„ of (2) into a geometrical 


Pa (cos), (K20), (70e) 


.P,, (cos $), (K20). (TOf) 
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series, just as in the plane case of §9. Thus we get 
instead of (62) for R, : 


ee CP (kya), OP (ka) on (ka) 
Ra=—} Pn(ka) * Palkia) ‘But? bn(k 2) 


æ L 
E(t Ry )REI+Ra,)- Í ne (71) 
As in § 9, we here again encounter the two factors 
(C-HR) Rall +R) = Ry.RERe, 
representing the attenuations caused by the first refrac- 


tion, the K reflexions, and the last refraction, while the 
factor 


Forai 
on (kaa) 
represents the phase retardation and attenuation as 
a result of the waves having travelled K+1 times to and 
fro between the centre and the surface of the sphere. 

It is of importance to note that the factors 


CP (haa)? (kaa) 

tend towards zero for both k,>0 and a->o, proving 
the fact that (a) for a totally reflecting sphere of finite 
size and (6) for a finitely conducting sphere but of 
infinite radius, and therefore tending to a plane, no 
reflected waves from the inside will be present. Limiting 
ourselves to the latter case (ao) we note that, even 
with the slightest absorption, [Jm(k,a)—>0], only I? , and 
O°, are left as components outside the sphere. The 
part O°, treated after the method of § 5, leads to the 
same result as obtained there, viz., the classical expression 
for the plane case in the form of the A integral for sec- 


12. Reduction of the Subsidiary Waves Of toa 
Sum of Residues. 


In §3 formula (7) it was shown, following Watson, 
that I, could be transformed into the sum of a closed 
contour integral and an open line integral along the 
imaginary axis. This inhomogeneous character of the 
result was due to the fact that the integrand was not 
an even function of n—}4. However, the splitting up 
of Io into an infinity of subsidiary waves O% (as affected 
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in § 11) enables us to write each of these waves as a closed 
contour integral only, encircling poles of the order K+ 2, 
so that the total field thus becomes rigorously a sum of 
residues, without any additional open line integral. 

This transformation is based on the following con- 
sideration. As each of the subsidiary waves Of (as 
given by (70d)) is represented by a sum over n, n being 
integral, and as we want to transform the sum into an 
integral with integrand, which, apart from the factor 
2n+1, is even in n—4, we are at liberty to continue 
the functions of n for non-integral values of this variable 
in different ways before we transform the sum into an 
integral. Incidentally the result thus obtained expresses 
each of the subsidiary waves O% itself again as an infinite 
sum of elementary waves Ok m where m denotes the 
number of times each of these waves has travelled 
round the sphere (compare § 3, equation (9)). 

In order to effect this plan we first notice the identity 


(EO (k a) E+. Pp (cos 3) = {2 (k aer) EH 
. P, {cos ($—Kr—r)}, 
with which, e. g., (T0 d) can be written 


¢ E (2n+1).g(n).P,{cos(3—Ka—m)}, (72) 
n=0 
with 
gin) = RRR. (e (k b)i Pk a) 


{nce ae ia (kr) 
C (ka) * OO (ka) 7 
As a consequence of the relation 
a) , a) 
C2) 19 (2%) =et™. nae (z) 

it follows that g(n—4) is built up of factors each even 
in n. This is also true for the “spherical reflexion 
coefficients ” such as Rs whereas the formerly used 
coefficient R (equation (2)) is not an even function of n. 

The fact that g(n—4) is even in n enables us, following 
the method of § 3, to transform the sum (72) first to the 
integral 
1 f ate ndn 


ES tot =i i Rarere 
“ ~œ +ie COS (N77) -9(m—})-Pn_12 {cos (3—Ka—2n)} 


i 
. (73) 
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As, further, along the whole path of integration Jm(n)>0, 
we may develop the factor {cos(nz)}-! as 


1 — 9 5 (—1)™. mmt, 
cos (N7) m=0 


so that (73) becomes 


w 
2 
Ok = a Ox, m? (74) 
m=0 
where 
— lje f> +ie s E 
k, m= ( : ) n dn gita(2m+t) -Ry2ReRo; 
t —2 +is 


: LOr 12( £15) 2 „o(ka) o 


EE aitaa E 
(Kaa) 
one ol kyr) 
. ~~~ .P,_3:2{cos (s—Km)}. . . (75) 
a (%10) ue { ( y 
This integral (75) can be closed by a half-circle of 
infinite radius above the real axis, so that its value can 
be represented by the sum of residues of poles n, of the 


order K+2. These poles are the zeros of M,_,=0, 
where 


My = [5 . Slog {2t(a)}| 


=k 


which occurs in (68 b) and is only slightly different from 
the expression N, defined by (6). Next making use of 
the asymptotic formula (12), each of the residues is seen 
to have as exponential part the factor 


n 2m} 
e (G+ 2m : 


showing, finally, that each of the terms Ok, » represents 
an elementary wave which has been reflected K times, 
refracted 2 times, and diffracted m times (having encircled 
the sphere m times) before arriving at the receiver. 


In considering the theory of the rainbow the most 
important elementary waves are O? , and O$ ,, whereas 
in the radio case the elementary waves I’ , , and O°, o form 
the essential part. A further analysis shows that the 
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sum I’, , +0%,,) can also be transformed to a sum of 
residues of poles of order unity, yielding 


dari = Ny. en 
It, o +02, o> (k ka)? ` nae (OM « i 
On Jn=ns 


e O 


and it is on this expression that we will later on base 
our numerical calculations in the radio case (see § 17). 


13. Transformation of the Elementary Wave O%, „ into 
a (2K -+ 8)-dimensional Exponential Integral, the 
Saddle Point of which determining a K-times internally 
reflected Goemetrical Optical Ray. 


In this paragraph we will develop a totally different 
aspect of our problem, leaving for the present the residue 
method. Here we will transform the general elementary 
wave Of m into a more dimensional exponential integral. 
It will be seen that the position of its saddle point deter- 
mines at once the only geometrical optical ray connecting 
the transmitter and receiver, which has been reflected 
internally K times. Moreover, a further development 
of the exponent round the saddle point provides us with 
an expression for the amplitudes of this ray. Performing 
this transformation for all values of K we thus succeed 
in writing the complete solution of our problem as the 
sum of an infinite number of geometrical optical rays. 

In order to affect this transformation of (75) we first 
note that the denominator only contains (@ and g 
functions, each of which possesses an exponential character 
which we can annul by multiplying the nominator and 
denominator each by 


Ca yolka) (on? yjolRea) y t, 


resulting in a denominator which only contains factors 
of the form ¢°) |. (2) . 01,2 (£), where the exponential 
character has disappeared. As the same was already 
true as regards the factors R,, Rẹ Ra, we are now 
prepared to transform each of the factors remaining 
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in the nominator into an exponential integral by the 
aid of the fundamental expressions 


1 OO ceca stay a 
Mied =o J earr- de, (18a) 
1 —r-ix i: 5 
nof?) pro ya : j a gents {TSD 
r-n æ% 


where —arg z<yn<7—arg z. 
Writing finally 


Pa -1:2 {cos (3 — Kz)} 
= i : i‘ en- 1/2) log 4 cos (9—Kr)}+isin (9~Er)}cosġ } dd 
Qa -F ? 


we obtain, instead of (75), 

"E (=1)" 

Em Sik 2a V br(2rika) t 

m g 
| _dnatm. f defar. -drxsee tm, (79) 
where the non-exponential factor A(n) is given by 
n. RRR 

pol byt) ED olka). ER aalkaa) oe aalkaa) H 
while the exponential factor exp. (fk, m) has as exponent 


A(z) = Z 


2K+4 
Si. m= ik ib cos 7,—ikya cos7,+1k,a 2 cos ty 
¥=3 


—th,@ COS T2K45 bth COS Tox 46 


2K+4 
tind Ti — Tat = TyY—ToK45 T ToKigt (2m+ 1e b 


+(n— $) log {cos (3—Kz)+¢ sin ($—Kz) cos ¢}. 
(80) 

Further, the paths of integration of Ti, 73... Tox; 
Tox+6 are as in (78 a), while those of 7, and 72x45 are as in 
(78 b). 

The physical character of (79) becomes at once apparent 
if we consider the position of the multi-dimensional 
saddle point of exp. (f m) which is given by 

af_ of af af _ af 


Ot, Or, `O rae On OD ”’ 
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the solutions of which are 
n=k,b sin 7,=-k,a sin 72=k,a sin 73= «.. =k SİN TK 4 


2K +4 
(Ta— T) + = ($r— Tte) t (T2K+5— Txx+—) = 3+2mz7, 
(81 5) 
MRS 228 ae Shs Fe ae A ee og it ve ey BS) 


If we now identify the variables + with the angles as shown 
in the schematic fig. 13, which we consider to lie in the plane 
¢—0, the expressions (81) determine exactly the geometrical- 
optical broken ray, which first leaves the transmitter at Q, 
is thereupon refracted at A,, and next internally reflected 
at A,,....Axy1, and which finally leaves the sphere as 
the refracted ray Axi2P. That this is the case follows 
immediately from the fact that (81 a) represents nothing 
but Snellius’s Jaws for each reflexion and refraction, 
and (81 5) expresses that the total angular deviation of 
the ray as measured at the centre O is the sum of the 
angles (T2— 71), ($T — Tə), etc. 

Having thus determined by the position of the saddle 
points the geometrical optical rays, the intensity of 
these rays follows from an approximate evaluation of 
the integrals (79), but before we proceed we first insert 
a short consideration on the physical meaning of the 
first, second, third, and fourth order approximations of 
integrals of the type of (79). 


14. The Physical Meaning of Approximations of 
Different Order of Saddle Point Integrals. 


A complex exponential integral, 
fed, 2 2 a n a (82) 


can often be approximated, as already considered, by 
first finding the position of the saddle point z,, determined 
by f’ (¢,)=0. The path of integration may then often 
be chosen in one of two ways, determined by 


(a) Im{f (z)}}=Imf{f(z)}, (steepest descend, Debye, 
Watson) ; 


(b) Ref{f(z)}—Re{f(z)}, (stationary phase, Kelvin). 
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In both cases the main contribution to the integral 
is furnished by the immediate neighbourhood of the 
saddle point z). In case (a) this follows from the fact 


Fig. 13. 


The geometrical optical ray joining Q and P. 


that the modulus of the integrand most rapidly decreases 
as we leave z,, while in case (6) the phase changes more 
and more rapidly as we leave 2), causing a cancellation 
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of the contribution to the integral by the parts of the 
path far away from z, *. 
Developing further in (82) f ane as 


z- ae ( zia 


Piai —f'"(z, + & HZD + T¥ (Zo) « 
[er dz =| dee i ard isa ohn 


we may get a first, second, third, and fourth approxi- 
mation to the integral by retaining of the exponential 
series one, two, three or four terms resp. only. 

Integrals of the type (82) often occur in diffraction 
problems. When they correspond to two-dimensional 
or three-dimensional space, the latter with axial 
symmetry, f(z) depends not only upon the complex 
variable z but also on two spatial coordinates x and y. 
In this case the first approximation diverges for all 
values of x and y when, as in our problem, the limits of 
integration extend towards infinity. However, this first 
approximation converges when the limits are finite, and 
is then usually called the Fraunhofer approximation. 

The second approximation to (82) leads to Fresnel 
integrals, whose physical interpretation is beautifully 
illustrated by the well-known Cornu-spirals. However, 
this second approximation may also diverge, viz., in the 
case when f” (z,)=0. This equation, considered as 
a function of x and y, determines a curve—the caustic, 
the envelope of the rays —where therefore this approxi- 
mation tends to infinity, as will further be shown in the 
next paragraph. 

Considering next the third approximation, which leads 
to Airy’s integrals, having the form 


| e-iz +iaz dz, 


we observe that they may diverge if simultaneously we 
have f(z) =f (z) =0, being two equations in v and y, 
and therefore determining isolated points in space only. 
These points are the cusps of the caustic, 2. e., the foci. 

Finally, the fourth approximation would only diverge 
if simultaneously we have f”(z) =f (2o) =f (z9)=0, 
which in physical problems would lead to three equations 
in the two variables x and y, which in general are contra- 
dictory. 

* For a further treatment of these questions we refer to the important 
paper by L. Brillouin, Ann, Ec. Norm. (2) xxxiii., Jan. 1916. See also 


Van der Corput, ‘ Compositio mathematica,’ i. pp. 15-38 (1934) and iii. 
pp. 328-372 (1936) ; J. Bijl, Dissertation, Groningen, 1937. 
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The multi-dimensional integrals (79) occurring in our 
problem are a more-dimensional extension of the integral 
of the type (82). Also in this case we can consider the 
first, second, third . . .approximationsas defined above, and 
they also diverge in a smaller region of the space coordin- 
ates the higher the order of approximation. Therefore 
a higher order approximation is necessary the nearer 
the caustic we wish to investigate the field. Especially 
in the case of the rainbow, for which at infinity the 
minimum of deviation coincides with the caustic, and as 
appreciable light is only to be found near this minimum 
we are forced to consider the third approximation 
(see § 16). Also in the radio case, between the trans- 
mitter and its horizon—at which latter point we reach 
the caustic—our second approximation remains finite only 
to diverge at the horizon itself. 

In the next paragraphs we shall consider in more 
detail the evaluation of the second and third approxi- 
mations for the integrals ofthe form (79). We saw already 
that they represent rays, and our next task therefore 
is to find their intensities. 


15. The Second-Order Approximation to the 
Intensities of the Rays. 
Our integrals of the form (79) being of a multi- 
dimensional type, we first give a short résumé of their 
second-order approximation. We therefore consider 


if] je efn 2 - -- tn) dzy, dzo . dza, (83) 
-0 Y- `- = 


the saddle point (z,°,2,°, . - - z,°) of which is determined by 


Now the second approximation of (83) is obtained with 
the aid of the more-dimensional Taylor development of 
f (% Z2- - - Za) round the saddle point, stopping at the 
second-order terms, leading to 


am (s) iral 
I~il, =en. z2, Meee zn) | | eee | 
— %0 


B Fp 
dro dz, day... dza. 
. . . (84) 
Phil. Mag. S. T. Vol. 24. No. 164. Suppl. Nov. 1937. 3I 


(zı —z1°}2 ( pa ) + (21- 22°}{z2—Z29} . { 
xe 
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Further, the exponent of (84) may be transformed to 
“normal coordinates ” with the aid of the linear trans- 
formation 


n 
Zi —2} = È pty, 
k=1 


n 
so that the exponent becomes 2 a; u,?, and the multiple 
k=] 


integral (84) reduces to the product of n complete error 
integrals. The result is 


iv)". | Pa | 


I, — ( . ef 22... zn?) 7 4 (85) 
Vaja. - sün 


where | p | means the determinant of the coefficients 
Pa- Moreover, according to a well-known theorem of 
linear transformations we have 


= A 
Vaaa, =| pal A/S, 


where A is Hesse’s determinant of order n of f(z), tken 
at the saddle point, and which is defined by 


ee 
a= zdz, 0 i 
Hence (85) leads to the final result 
. Qn)" 
I,= (iV 20) Bf Z 22% > -+ Z£) 
2 V/A , 


where, as we shall see further on, the first factor gives 
the amplitude and the last the phase in our ray problem. 


(86) 


We are now ready to apply (86) to (79), in order to 
find the second approximation of (79). We therefore 
identify f(z,, Z2 - - - Za) of (86) with fk, m as given by (80), 
and this for the values of the z’s in the saddle point which 
are known from (81). Referring to fig. 13 we can express 
the result as 


(fix, mlo=t{k, QA, +h (AA, +. -. + Ax Ay) 
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< 


where Sx is the “ eiconal ” * belonging to the ray which 
has been reflected internally K times, and which expresses 
the phase of the ray. In order to find its amplitude 
we must evaluate 4, and a further consideration shows 
that 4 is most elegantly expressed with the aid of 


es stk T 1 2K +4 l 
ðn — ke, bcost, | kyacost, v-3 k,acos Ty 
1 1 
ae (88 
as 
A=(—1)*.k,b. kya. (ka) =+. ka. kyr. (n—4) sin $e- 


os 


- COST, COS To. . . COST. <a e 
1 2 2K+6 on 


The geometrical meaning of a is obtained if, by keeping 


On 
the transmitter Q (fig. 13) fixed, we imagine the receiver P 
to move a little distance along its circle round O. Thereby 


3 and n vary with a differential quotient ad . It is 


further of interest to note here that the physical geo- 
metrical meaning of the variable n, which originally 
appeared as a separation constant in the harmonic 
series (1) from which our whole investigation started, 
according to (81) is nothing but 27 times the length of 
the orthogonal from the centre O on all cords of the broken 
ray (expressed in wave-lengths), and which is even true 
for those parts of the ray which lie outside the sphere. 

Returning to (79) A(n) has now to be taken at the 
saddle point, and further introducing in A(n) approxi- 
mations of the type 

1 
(1) C) oae ee fe. a 

Ka sin T° ON sin 7 (ka)? COS T (89) 
and limiting ourselves to m=0, we find for OẸ o the final 
result (second-order approximation) 


5 i. RaRa R oa y V kya sin Ta ei8k (90) 
ES / . a OF VWcosT,.CcOSTy., 
kb kyr sin 3. On -K+ 


* See also. Ch. Manneback, ‘ Recueil des travaux de l‘Assemblée 
Générale de |‘ Union Radio Scientifique Internationale, tenue à Londres 
en Sept. 1934”: Bruxelles, 1935. 
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This formula, which represents the ray of order K, 
confirms the considerations given in § 14 in so far as it 
becomes infinite where a =0, i. e., on the caustic. 

With a view to the radio application, where the ray 
of order —1, which is the ray which has been reflected 
only once externally and which has to be computed 
separately, we give here its approximate expression in the 
form 


vka sin To Ry 5-1 


Oio ~ a8 
ikyb.kyr \/ sins. By -COS Ta- C08 4 


e*1(Rit+Re) 


= R pra neem 
Maa 1° ER, ARa) 


(91) 
where 
Cr Rovia sin 79 
ae FA E 


which, after a slight transformation, becomes 


N a(R +Ra) V sin 72.008 T3 
11°” Vbrsin $(Ryr cos r,+Rbcos 7) 


Summarizing, we refer to fig. 14 for the meaning of 
the variables 7 and R. 

Formula (91) tells us the interesting physical fact 
that the reflected ray is not given by the simple product 
of the primary ray and the reflexion coefficient R,,, but 
that a factor «,,, which we call the “ divergence factor,” 
also enters. This is physically obvious, because the 
reflexion in our case occurs against a curved surface, 
so that the divergence of two nearby primary rays, 
contrary to what happens for a reflexion against a plane 
boundary, is increased, and therefore the amplitude 
of the secondary ray is diminished, due to the increased 
divergence (astigmatism). This is confirmed in the two 
special cases (a) r=a (receiver on the sphere), where this 
extra divergence obviously disappears, and we therefore 
obtain «,,—1, and (b) aco (reflexion against plane), 
where, again, «=l. 


(92) 
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Applying, finally, our second approximation to I, , 
(of (70)), we find 
oR 
I, oo r= ik R’ 


Fig. 14. 


Q 


The direct ray together with the once externally 
reflected ray. 


so that in the radio case the second approximation to 
the total field becomes 


eit R 
Tgp ~T 2,0 +081, o~ TRU +a R pet, (93) 
1 


This extended reflexion formula differs from one occasion- 
ally found in the technical literature, where the divergence 
factor «,, is taken as unity, and instead of our spherical 
reflexion coefficient R,, the ordinary plane reflexion 
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coefficient is used. However, the numerical difference 
is small in most radio applications. 


16. The Third-Order Approximation to the 
Intensities of the Rays. 


Having in the previous paragraph considered the 
second approximation of the elementary waves Oj,» as 
given by (79), we shall now investigate the third approxi- 
mation of these waves. As pointed out already in §14, 
the third approximation is the fundamental one in the 
theory of the rainbow, because here light is practically 
only present near the minimum of deviation, which 
coincides with the caustic at infinity. 

Developing the exponent (80) round the saddle point 
up to the third-order terms, it is found that the variables 
r occur separated, so that the multiple integration over 
all the rs can be affected first, thus leading, with 
tTy—Ty°=u, to a product of integrals of the type 

POD pede A 
i) g7 0 Fiu? iyu dy — 6 Zi at” 3a 


ae 
3V a 


3324 a faA* 
where the paths of integration are as given in fig. 15, 
W, belonging to the functions (2) and W, to the functions 
(2), 

The last approximation of (94) can be obtained in two 
different ways: (a) as the first term of the asymptotic 
series of the Hankel function, and (b) by applying (86) 
to the integral of (94). In both cases it is assumed 
that 


>> a, 
8> >ay, 
which in the case of (79) comes to 


Having thus affected all the 7 integrations we are left 
with the integrals over ¢ and n. The integral over ¢ 
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appears to be a nearly complete Fresnel integral, which 
can be approximated by unity. The remaining integral 
over n contains an exponent which can be developed 
again up to the third-order terms in »—%p , thus leading to 
an integral of the form (ef. fig. 15) : 


6 W fe 2 o 
| a a gv oe ER 33/2 ie) ; 
Z wa 


Path of integration of third-order integrals. 


giving as final result 


Š ir+iS 
k o p ary Gain Scos r. coame RaRa 
In -2(89)9/(28)? ra! 
ge Te Oe eel 


an? 1 On? 
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where in the result analogous to (88) we have, apart 


3 ; 2 . : 
from a , also introduced , the latter being given 
by the geometrical relation 
as sint, sin Ta 2K+4 = gin Ty 
Ən? k B® cos*z, "` ky?a*cos? t, vos ka? cos? ty 

sin T2K 45 sin T2K +6 


kžažcos*rmgs krt cos ryo OO 
(95) therefore constitutes our third approximation to 
the elementary waves represented by (79). 


The factor Hf} contains in the argument the factor 
3 
pa , which is zero on the caustic and has different 


signs on its sides, and this nicely demonstrates a well- 
known general property of waves near caustics, viz., 
that their amplitude falls off exponentially at one side 
of the caustic and in an oscillatory manner at the other 
side, this fact being demonstrated by the exponential 
behaviour of H®3(x) for v>0 and the oscillator character 
of the same function for <0. This behaviour is due 
to the fact that at the convex side of the caustic, where 
as 


On <0, each point is in general reached by two different 


rays causing an interference pattern. 

It is also possible, as a further analysis shows, to give 
the geometrical meaning of the argument of H, in (95). 
In fig. 16 the caustic is represented by the full-drawn 
curve C, while the wave-front W, cutting the caustic C 
in the cusp W, orthogonally, is shown dotted*. Further, 
there are represented two rays, W,P and W,P, touching 
the caustic at K, and K, resp., and which rays arrive 
at P with a path difference A—W,P—W,P. Due to 
the caustic being the evolute of the wave-front W, this 
expression may be replaced by 


4—K,P+PK,—arc K,K,, 


* See also R. Wood, ‘ Physical Opties,’ p. 148 (New York, 1919). 
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and a somewhat elaborate transformation, moreover, 


shows that 
323 
a (3n) 


“seme 
an) 


, 


Fig. 16. 


Caustic with wave-front. 


so that the argument of H; in (95) simply becomes 
4.e€*7.k, ^A. The asymptotic expression of H., then 
yields two terms, representing the two rays with a phase 
difference k. We will return to these third-order 
approximations in § 18. 
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17. Numerical Results for the Radio Case. 


Summarizing, we have now available for numerical 
interpretation for the radio case the following expressions :— 


(a) Watson’s expression, which is a special case of (15), 
viz. for b=a, and therefore only refers to the case where 
both transmitter and receiver are situated on the ground, 
and which is based on the residue method. 


(6) The expression (15), which is also based on the 
residue method, but here either the transmitter or the 
receiver is situated above the ground *. 


(c) The expression (77), covering all cases, and which 
was derived with the residue method after the total 
field was broken up into elementary waves and after 
it was shown that in the radio case waves reflected against 
the centre of the sphere were negligible. Thus the 


denominator © of (77) differs from Watson’s a of (15) 


as a consequence of the fact that the direct radiation 
and the only once externally reflected radiation were 
maintained. 


(d) The reflexion formula (93), which was derived 
after the total field was broken up into elementary 
waves, the direct radiation and the only once externally 
reflected radiation being maintained, the derivation being 
obtained with the aid of the saddle point method, 
yielding (1) the spherical reflexion coefficient R,,, and 
(2) the divergence factor «,,. 


As to the ranges of applicability of the different 
expressions, these are most clearly described when we 
consider first the field far beyond the horizon of the 
transmitter. There formula (77) can be used, and its 
first term already gives a very good approximation. 
When we approach the horizon from beyond the number 
of terms of (77) to be considered rapidly increases, 
although at the horizon no more than six to ten terms are 
necessary in practice. Going further towards the trans- 
mitter, as we have passed the horizon, the number of 
terms of (77) to be retained rapidly increases still further, 
and in our examples we have gone up to eighteen terms. 

* After the present investigation was finished a paper appeared by 


Wwedensky, ‘ Technical Physics of U.S.S.R.’, iii. p. 915 (1936), giving 
a formula analogous to our expression (15). 
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This shows the practical applicability of our general 
expression (77), which, however, is formally valid every- 
where even when the receiver approaches the transmitter 
indefinitely. 

We have further found it convenient in applying (77) 
to limit ourselves to two limiting cases (which have 
already been considered before), viz., 


(1) Total reflexion : c= (see §75), 
(2) Maximum absorption : 
a kè] k? 

[è| = s113 ° ./p 22 1 

H T yka j k 2—1 
(see § 8), 

for in these cases the values of 7, and 7,’ can easily be 
calculated (see (40 a) and (55)). 

It is thus seen that, while approaching the transmitter, 
the number of necessary terms of (77) soon becomes 
unwieldy. Mostly there is still left a transition region, 
dotted in the figures to be given later on, but then we 
soon approach shorter distances from the transmitter, 
where we can apply the reflexion formula (93). The 
values of the field in the transition region can, however, 
mostly be properly estimated by interpolation. 


<<1 


Already in Part I. of this paper we gave some numerical 
results also for either the transmitter or the receiver 
above the ground, these cases being equivalent according 
to the principle of reciprocity (see fig. 8). The effect 
in (77) of raising the height is expressed by the factors 


ge {a 

a 0 (kya) n=z+a Agg 
present in each term of the series. If both transmitter 
and receiver are raised above the ground two of these 
factors are present, the second one with b replaced by r. 
In order to discuss the numerical effect of these factors 
f, the ¢ functions can be replaced by the approximation 
(17), yielding (after a slight simplification based on 
z'83,5<<1, where $,=cos1b/a) : 
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and a similar expression when 7, is replaced by 7,’. Here 
we have introduced the new variable y=x'* 3=(ka)"8; 
the special value of y for a point at the horizon is 
designated by y,=x"%3, if either the transmitter or the 
receiver is above the ground. If both are raised above 
the ground (see fig. 17) the transmitter as well as the 
receiver each has its own horizon, for which the y’s 
are denoted by xy,;=2'?S, and x2.=2'*S,. 


Fig. 17. 


Definitions of 3, and 3». 


Returning to (97) a very important physical property 
is expressed by these factors, viz., the increase of the field 
with the height. This increase appears to be a function 
of (a) the numerical constants 7, depending upon the 
electrical constants of the ground and the wave-length 
and (b) the geometrical variable y 9. Now it appears 
that the argument of the Hankel function in the nomin- 
ator of (97) is of the order 1/3 x,3, and the numerical 
behaviour of this Hankel function is totally different in 
the two cases 1/3 y,?<1 and 1/3 x,°>1, corresponding 
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resp. to h< 56)? and h>56A?%, where h=b—a and A are 
both measured in metres *. 

When h<<56A?'* the Hankel functions in (97) may be 
approximated by the first terms of their power series 
development, leading to 


fa! 4-0-459eBi (y 2—27) —0-0417 (x —2r;}?. .. 
* ~ “T+0-459e2**(— 27,)—0-0417(x92—27,)%.-. 

When, on the other hand, h>>56A?/, we approximate the 

Hankel function in the nominator by the first term of 


the asymptotic development, while the Hankel function 
in the denominator must first be written as 


(98) 


Hz) = 


m {J lz) +3 alih, 


and then the first term of the asymptotic development 
of both J,,, and J_,,, can be taken. This leads to the 
following final formule, valid therefore for x°>>3, 
eas x2>>3: 


Values of 7: 


(1) o= 0: T,~4{3n(s+4)}°* i em 
(2) [8] <<1: rAR{3a(s +3). c=? 8, 


( for s<< z) i 
Values of height-gain factors : 
(1) o= 
(— D- 275) ft (= E 


m~a ay 


. 2V Xo 


(2) [8] <<1: . (99) 
k (1 3 2r,’ Brat T 7 
12 (—1)'(—27,')"* | 


hae = “Ti ex 9 


* In an investigation of the same problem by Mr. T. L. Eckersley. 
who recently kindly communicated to us his results (not yet published) 
obtained with his “ phase integral method,” the effect of raising the 
transmitter is called the “ height gain ” ; he also finds a critical height 
given as h=50A?/°, which corresponds well with our A= 5642/3, as found 
above. 
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(100) 


en 


Attenuation factors for b=r=4 : 


~V rix. 2 — 
et o 282, 4/2 wig. se a| 


tot 


(l) c=oo: 
pr 

(2) [8;<<1: >~ 
[3|<<1: 7 
Attenuation factors for b4a, r=a 


(1) c= 
Th E SODH rart- 
o | (—2r,)**° 
(101) 


a peer 
Hio 24 2 ary 1 9 T4 
Few EV i 
xè a at, 3/2 

Xo? 


1T,'x+ 


h> 
et 
© 
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D 
ag 
© 
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3 
a 
o 
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œ 
th 
Q 
Q 
| 
bee | 
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(1) e=: 


TT eot TX 
~ iXX, 
274 \3/2 + Be — 5)" 


if. 
ix? 
aT sx 3 (1 We 
? 


€e 
2 f V —2r, 


(2) [è]<<1: 


Io 346 82 7 / TX 
vent EG NE (1 ea 


(102) 


(—2r; y 2.e : 
> {i 96 5 ra V? | 
i 
Figs. 18-20 were E with the aid of the above 
formule, where occasionally a still higher approximation 
. These figures 


for f, than given above was made use of. 
h,=0, 
h,=0, 


refer to 
Fig. 18, o=o, h,=100m., 
Fig. 19, o=10-8, «=4, h,=100m., 
o=10-%, «=4, h,=h,=100m., 


Fig. 20, 
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Influence of the wave-length on the field directly on the surface of the 
earth for a transmitter at a height of 100 m. (s=o). 
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Influence of the wave-length on the field directly on the surface 
of the earth for a transmitter at a height of 100 m. («e=4; o=10-). 
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and they represent | JZ, | divided by the modulus of 
the primary field e“:®/(¢k,R), and therefore represent 
the attenuation factors with which the undisturbed 
radiation must be multiplied due to the presence of 
the earth. 

Fig. 18 (for =o) clearly shows that waves longer 
than about A=10 em. do not present a pronounced 
effect at the horizon, while shorter waves rapidly fall 


Fig. 20. 
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Influence of the seis on the — ei on transmitter 
and receiver at a height of 100 m. (e=4; o=107%*). 


off behind the horizon, tending to a discontinuous drop 
to zero for waves A>0. 

Fig. 19 (with absorption) shows a similar behaviour 
of the waves, apart then from the absorption which is 
already present for all wave-lengths considered. The case 
of the limiting wave-length A-0 is here fully determined 
by «e only, the conduction current becoming negligible 
in comparison with the displacement current. 

Fig. 20 (with absorption) again shows a similar behaviour 
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of the waves ; only the wave A=70 cm. has a pronounced 
minimum before the line of sight. The interference 
maxima and minima, which undoubtedly are present, 
fall outside the part of the curves so far calculated. 

Herewith we conclude our investigation of the radio 
problem, and in the next paragraph we consider in further 
detail the application of our formule to the theory of 
the rainbow. 


18. Theory of the Rainbow. 


While in the radio case it was sufficient to consider 
the direct wave together with the wave which was 
once externally reflected only, because of the big absorp- 
tion of all the other waves entering the earth, in con- 
sidering the rainbow the conditions are just reversed, 
because the rainbow is observed near the minimum of 
deviation of rays which have been once (K=1; primary 
bow) and twice (K=2; secondary bow) reflected inside 
the raindrops. 

But as our theory has been developed generally, we are 
ready to apply it in the present case. It only needs 
two extensions: (a) the source (here the sun) will be 
considered at infinite distance, and (6) the receiver (the 
observer) also at a very great distance. Moreover, if in 
the radio case we consider a vertical electric dipole at 
an infinite height above the earth, due to the radiation 
characteristic of this dipole being zero in its own axis, 
no energy would reach the earth at all. We now have 
therefore, following Debye*, to replace the electric dipole 
of the radio case by a combination of a contribution 
of an electric and a magnetic character, the sum of 
which represents a plane polarized wave. 

Whereas in the classical theory of the rainbow the 
usual procedure was first to find the wave-front after 
K internal reflexions, then approximating this wave-front 
round its inflexion point and summing the different 
contributions from there with the aid of a Kirchhoff 
integral, leading in its turn to Airy’s theory (the so-called 
Rainbow-integral), we will here follow a totally different 
way, and develop the result starting from our rigorous 


* P. Debye, “ Der Lichtdruck auf Kugeln von beliebigem Material ”’ 
(Inaugural Dissertation, München, 1908); Ann. d. Phys. xxx. p. 57 
(1909). 


Phil. Mag. 8. T. Vol. 24. No. 164. Suppl. Nov. 1937. 3 K 
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solution of the wave equation, which was split up in the 
subsidiary waves O{ of (70d). The third-order approxi- 
mation of the multi-dimensional saddle point integral of 
Ox, as considered in § 16, then leads to an extension 
of Airy’s classical results. It is just this extension 
which here enables us to determine in a rigorous way 
(a) the absolute intensity of the rainbow, (b) the losses 
at the reflexions and refractions, and (c) the state of 
polarization. Whereas these points had formally to be 
introduced ad hoc in a more or less arbitrary way *, 
they are quite naturally present in our theory, due to the 
fact that it was derived from a rigorous wave solution. 

In order to be able to compare later on our results 
with the classical ones we shall now recall some formule 
of the older theory, which gives us the opportunity to 
compare the notations used here with the classical ones. 

Referring to fig. 21, let the sun be situated at the 
infinite point Q, corresponding to 9—0, and the observer 
at the infinite point P at an angular distance 3 from Q. 
A typical ray belonging to the primary bow (K=1) is 
represented by QA,A,A,P, this ray having a total 
deviation D=az—S$. As before, we will designate the 
rays by the value of the parameter n, see (81 a), which 
is given by 
n=k a sin ra= ka sin T3=. . . =k SIN 72K 44=h asin Tak 45- 


Expressed in the variable n, the deviation D(n) then 
becomes 


D(n)=2 sin-( 7) —2(K +1).sin-1/ ig) +, (103) 


and the minimum of deviation, round which the light 
energy is concentrated, is given by 
dD(n) 
dn 


Further, for this ray, suffering the minimum of deviation, 


it is easily found that the corresponding 7, becomes 
T2 min determined by 


=0. 


_ fp 2 
SiN 72 min = ae eed 7 La - « (104) 


K({K-+ 2) 


* See, e. g, W. Möbius, Ann. d. Phys. xxxiii. p. 1493 (1910), and 
xl. p. 736 (1913). 
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from which it follows that the ray of lowest order corre- 
sponds to K=1 (causing the primary bow), because for 
water in the optical range we have k,?/k,2=1-78. The 
secondary bow corresponds with K=2. 

Referring further to fig. 21, the classical theory 
considers the wave-front W determined by 


and the ray of minimum deviation cuts this wave-front at 


jts inflexion point O’. It was further usual to introduce 
new rectangular coordinates € and y with O’ as centre. 


The geometrical-optical ray QP in the primary rainbow, 
together with a wave-front W. 


Expressed in these coordinates the first approximation 
for this wave-front round its inflexion point was given as 
— h 3 

a 


with 


p= KAE +2)? V(K+1)?—k,?/k,? (105) 
(K+? ` (kjk P1 ? 

and it was the interference of the waves originating 

from this wave-front round O’ which lead to the classical 

theory. 


3K 2 
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In the present electromagnetic theory of the rainbow 
we have to start with a plane polarized wave as primary 
field, 7. e. 

'k 2c 


E,=—tk,.e-**; H,= PE. E, =E, =H, =H, =0. 


(106) 


This plane wave is polarized in the y direction, but 
we shall later on have to take the mean square values 
over all directions of polarizations at right angles to 
the z-axis in order to account for the unpolarized light 
from the sun. The two parts of the Hertzian vector, 
introduced by Debye, necessary to describe this field are 


Trot = Tot, el + Tot, magn? 


and the components E and H of the electromagnetic 
field are derived from this wave potential as 
, y 


E=curl curl (ra) + a . curl (rIImagn), | 

=: okt ( (107) 
= Sua curl (ra) +curl curl (rI magn)- l 

The proper primary field (106) is then obtained from 


(107) by taking pr =Hpr, a H pr, magn, Where 


aee 
I ac ikr sin 3 | 
l (e renion (kır)+isin (k,r) cos 3}, . (108) 
a _ ic sin $ 
pr, magn ™ pw ; rsin j 


. {e7 #1" 8 3 __ cos (kyr) +7 sin (kyr) cos 3}. 


Through this artifice the Hertzian vector has a radial 
component only given by tIio- 

In order again to fulfil the boundary conditions we 
want instead of (108) the equivalent harmonic series 


Th ve: e1. = — COs $ 
. E(2n+1)i-Yq (ker) . By (cos $), 
n=l 
Il, magn ~~ —sin $ i 


i 5 (ant 1)i-*,(kyr) . P7 (cos 3). 
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As is also well known in the plane problem the 
boundary conditions are different for Ia and Imago- 
So far in this paper we have used the electric boundary 
conditions only, whereas now we have to fulfil the two 
separate conditions 


s (rIlhot,ea) and kIT... continuous at r=a, 


0 f 
ar (72D sot, magn) CD. Liot, magn CONtinUOUS at r=a. 


Extending the results of § 10, where we derived the 
spherical reflexion coefficients R, and Rẹ which we 
shall further call Ria and Rze, we shall have to 
introduce as well Rii, magn and Rye, magn, Which are at once 
derived from the magnetic boundary conditions, and which 
are found to be 


Rn, magn — 


~[agmeermen] + [egotoe rote] 


z=kga 
Ry», magn = l+ Ryy, magn. 

All the results obtained heretofore for the electric 
dipole may at once be applied in slightly extended form 
for the present rainbow case. The essential differences 
are : 

1. Here we have Ha-+Hmagm instead of Ma only. 

2. Both the source Q and the observer P are 
at infinite distance from the centre of the sphere 
(raindrop). 

3. The occurrence of $$}. Pz! (cos $) instead of 
P? (cos $). 


Now the rainbows of different order K are rigorously 
given by the subsidiary waves O¢. Thus we find instead 
of (70 d) the series 


(0%) a= —4 cos ¢. Z (2n+1)i—*. (Ri2R22Ro1)e1 


EDk) EH Nka i 
l E : o os i CO (kr) : P- (cos 3), 


and an analogous series for the magnetic contribution. 
As, moreover, the actual bows are concentrated near 
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the minimum of deviation (i. e., the caustic), we have, 
as explained in § 16, for numerical evaluation to use at 
least the third-order approximation for the elementary 
waves Of ,, for which the electric part is easily found to be 

Sx a 
| Oc. Pein 
\ x, odel tk, . 


FE = i K 
6k,a.sin7,.sin 3' (Ri2R22R21)ai . COS $ 


- D 2:/dD\3 / fd2D\2 Fa i 
dn la)/a) oa] Sin) a 

T - this aye 3 (110) 
dn? 3 | an 


which is the equivalent of (95) and where the derivatives 
of D(n) are at once found from its definition (103). The 
expression for the corresponding magnetic part is again 
quite analogous. 

In order, further, to obtain the actual field components 
from the wave potential (110) we apply (107), bearing 
in mind that only Sg has to be differentiated, because 
ka>>1. Now the physiological intensity of the observed 
light is usually considered as given by the E vector * 
and can therefore be represented by I=E2, where the 
dash over the E? means the value averaged over all polar- 
ization directions 6. Now E,=0, and we are left with Eg, 
corresponding to the component at right angles (1) to the 
are of the rainbow, and Eè, which corresponds to the 
component along the arc (it) of the bow. Calling further 
I, the intensity of the E, vector of the incident light, 


we are interested in 
I _ ES e= 
D Ex Hyp 
I, _ Es? us) 
I, Ej?\” H 
We thus obtain from (110): 
(a) the absolute intensity 


i, 7 (Ry2R22Roy).? 
1, 12° kya.(k,r)?.sin Ta sin’ D 


dD\? ( dD\? 2 

(PV fal 2) 

2 mes (eb » AN) 
a Ga) 


* M. Born, ‘Optik’ (Berlin, 1933), p. 24. 


